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Abstract - Two-dimensional orbital compass model is studied as an interacting itinerant electron 
model. A Hubbard-type tight-binding model, from which the orbital compass model is derived in 
the strong coupling limit, is identified. This model is analyzed by the random-phase approximation 
(RPA) and the self-consistent RPA methods from the weak coupling. Anisotropy for the orbital 
fluctuation in the momentum space is qualitatively changed by the on-site Coulomb interaction. 
This result is explained by the fact that the dominant fluctuation is changed from the intra-band 
nesting to the inter-band one by increasing the interaction. 



Introduction. — Orbital degree of freedom is one of 
tjhe central and not yet solved issues in correlated electron 
£H systems PQE]. This degree of freedom is recognized as 
O a key to elucidate several exotic phenomena in solid state 
, ^ p hysics such as colossal magneto-resistance effect [3J , iron- 
Based high-T c superconductivity [3], and so on. Essence 
i-Hof the orbital is a directional nature; electron motions, 
^ and electron-electron interactions depend on directions in 
a crystal lattice [5] . This characteristic brings about sev- 
eral non-trivial results such as a macroscopic number of 
degeneracy in the orbital configurations, the dimensional 
^ reduction in the effective interaction and others [SHE] • 

, One of the well studied orbital models is the Kugel- 
Khomskii model [3] , which describes the interaction be- 
. . tween the localized ions with spin and orbital degrees of 
freedom in the nearest neighbor (NN) sites. The inter- 
^ actions on each NN bond are represented by products of 
?_ithe Heisenbcrg-type spin part and the orbital part, which 
5^ip described by the pseudo-spin (PS) operators, I*. In 
the latter, the orbital interactions explicitly depend on di- 
rections of the bonds in a crystal lattice. This model is 
obtained from the Hubbard-type itinerant electron model 
with orbital degenracy in the strong coupling limit. It 
is recognized that the magnetic and elastic properties 
in several transition-metal compounds, e.g KCUF3 and 
LaMnOa, are well reproduced by this model Hamiltonian. 

Another well known and simple orbital model is the or- 
bital compass model, where the orbital degree of freedom 
is only taken into account, instead of the spin-orbital en- 
tanglement in the Kugel-Khomskii model. In this Hamil- 



tonian, the component of the PS operator concerned in an 
interaction explicitly depends on the direction of a bond, 
e.g., only the x component Tf for the interaction along the 
x direction. In particular, the compass model in the two- 
dimensional square lattice has been studied so far inten- 
sively and extensively from the view points of the orbital 
order in a Mott insulator [TO] as well as qubits in the quan- 
tum computer [TT], and the topological order [T2]. Some 
new concepts, such as the directional order |13j . a non- 
trivial dilution effect [2], and the generalized Elitzur's 
theorem [T5] , have been proposed through the theoretical 
examinations. 

In this Letter, we study the two-dimensional orbital 
compass model from the view point of the interacting 
itinerant electron system. We identify the Hubbard-type 
model from which the orbital compass model is derived in 
the strong coupling limit. The model Hamiltonian is ana- 
lyzed by the random-phase approximation (RPA) method 
and the self-consistent (SC) RPA method. It is found that 
anisotropy in the orbital fluctuation strongly depends on 
the on-site Coulomb interaction. The results are inter- 
preted by the intra-band and inter-band nestings. Rela- 
tions between the present results and those in the original 
compass model are discussed. 

Model. — We start from the orbital compass model 
on a square lattice in the x-z plane defined by 
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where J is the exchange constant and < ij >i represents 
the NN sites along the /(= x, z) direction. The PS op- 
erator Ti with a magnitude of 1/2 describes a doubly- 
degenerate orbital degree of freedom. This is given by 
Tj = | 5Z 7 y c| 7 <r 77 /Ci 7 ', where c^ 7 is the annihilation op- 
erator for a spin-less fcrmion with orbital j(= a, b) at site 
i, and cr are the Pauli matrices. 

Here we introduce the Hubbard-type spin-less fermion 
model from which the orbital compass model in eq. (QJ 
is derived by the pcrturbational procedure in the strong 
coupling limit. This is given by 
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Fig. 1: (a) The band structure in the first Brillouin zone. 
The two bands contact with each other at four points of k — 
(±7r/2, ±7r/2) and (±7r/2, =F7r/2). (b) Fermi surfaces in the 
k x -k z plane. The red and blue regions are the electron- and 
hole-pockets, respectively. Black and white arrows represent 
the inter-band nesting at q — (0, it) and the intra-band nesting 
at q = (-7T, 0) respectively. 
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where rii 7 (= 
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cj 7 c; 7 ) is the number operator. In the first 



term in eq. ([2]), tj 1 is the transfer integral between the 7 
and 7' orbitals along the direction I, and /i is the chemical 
potential. In the second term, U is the on-site Coulomb 
interaction between the spin-less fermions with different 
orbitals. The matrix elements of the transfer integral are 
determined to reproduce the compass model in eq. ([I]) by 
the second-order perturbation with respect to the transfer 
integral. First we consider the exchange interaction along 
the z direction, i.e. T?T?. This Ising-type interaction 
requires a condition that the transfer integral along this 
direction is finite only for one of the two orbitals, i.e. 



t z = \t{l + a*), 



(5) 



where t is a constant. In the same way, the transfer inte- 
gral along the x direction is given as 



(6) 



We consider the so-called half-filled case where a number 
of fermions is equal to that of the lattice sites. 

The transfer integrals introduced above give the follow- 
ing energy bands in the momentum space, 
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with 



Sk = —t(l + <J X ) cos k x — i(l + a z ) cosfc z 



(7) 
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We take a lattice constant as a unit of length. It is worth 
to note that these band-dispersions are generalized as 

£fc(0) =-*[! + er(0)] cosk x - i(l + a z ) cos k z , (9) 



where a (6) = a z cos 9 + a x sin 6 with a real number 9. It is 
obvious that the bands with 6 = tt/2 are reduced to those 
in the compass model. The parameter values 9 = 2ir/3 and 
7T correspond to the models where the doubly degenerate 
(3z 2 — r 2 jx 2 — y 2 ) orbitals and the (yz/xy) ones are intro- 
duced at each site in a square lattice, respectively [TBlf!?] . 
By diagonalizing the matrix, we obtain the two bands 



=(±)fe 



(10) 



where = — t(cos k x + cos k z ) ± ty/ cos 2 k x + cos 2 k z are 
the band energies and dk n is an operator derived from Cfc 7 
by the unitary transformation. These band dispersions are 
shown in fig. [TJa). The two bands touch with each other 
at the four points of k = (±7r/2, ±7t/2) and (±7r/2, =fw/2). 
There is the particle-hole symmetry in the half-filled case. 
As shown in fig. [TJb), the system is a semi-metal, and 
both the electron- and hole-Fermi surfaces are squares. 
The square-Fermi surfaces are maintained in the bands 
with an arbitrary value of 6 in eq. ©■ Perfect nestings 
occur at q = (n, 0), (0, it) and (it, it). 

Method. — The model Hamiltonian in eq. (|2j with 
eqs. © and (gj is analyzed by using RPA and SC-RPA. 
We introduce the one-particle Green's functions defined 

by 



r/3 

G^(k,ie n ) = - dr(c fc7 (r)4 y )e 4£ " T , (11) 
Jo 

where e n = (2n + 1)itT is the Matsubara frequency, r is 
the imaginary time and /3 is the inverse temperature. The 
Green's functions at U = 0, termed G°(k,ie n ), are given 
as 
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Fig. 2: Diagrammatic representations of (a) the self energy 
£(k,ie n ) and (b) the effective interaction V(k,k';q). The in- 
teraction U is defined in eq. (|20p . The dotted lines are for the 
one-particle Green's functions as the external lines, and the 
bold lines are for the full one-particle Green's functions. 



with coefficients 



A k = ( 1 T cos k z /yj cos 2 k x + cos 2 k z ) /2, (15) 



and 



Bk = — cosk x / \ 2\J cos 2 k x + cos 2 . (16) 
We also introduce the orbital susceptibility defined as 
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) and?^ 77 } = Efc4 7 c fc+<ry'- 
From now on, for simplicity, we use the abbreviation for 
the suffix as a = (7,7') which takes f= { a , a), |= (6, 6), 
+ = (a, 6) and — = (6, a). Furthermore, we define the (zz) 
component of the susceptibility as 
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In RPA, the susceptibility is given by 
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and x° is the 4x4 matrix of Xa/3(<Zj * w n) which is the bare 
susceptibility at [7 = 0. 
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Fig. 3: Contour plots of the (zz) component of the orbital 
susceptibilities x5f A as a unit of t -1 obtained by the RPA 
method. Temperature is T/t = T c RPA /t + 0.01. Parameter 
values of U/t are chosen to be (a)0, (b)2, (c)3, and (d)5. 



In order to consider the band-modification effects, be- 
yond RPA, the susceptibility is calculated by the SC-RPA 
method, where the one-particle Green's functions and the 
susceptibilities arc calculated self-consistently. The sus- 
ceptibility is given by the same form with eq. (|19p where 
X° is replaced by x° defined by 



X(7i72)(7i72)( 9 ' ia; ™) 



2^72 7(( fc ■ 



q, ie m + iio n )Q l2ll (fc, ie m ). 

(21) 



The one-particle Green's functions are obtained by the 
Dyson's equation Q = [(Q )^ 1 — S]^ 1 where the self-energy 
is given by 

£ 7i7 ,(fc,i£n) = T Y1 I ?Sf2 y (7l 7^(7172) (9-^™) 
m •'B.Z. V ) 

x g 72ll (k + q,ie n + iuj m ), (22) 

as shown in fig.[5Ja). Here we define the effective interac- 
tion V{q) by 



v = (i-u Xo y 



l U--Ux U. 



(23) 



Schematic diagrams are shown in fig.f2Jb). The last term 
in eq. (|2"5|) is required to avoid the double counting of the 
diagrams. We note that, because of the on-site interaction, 
V does not depend on both the momenta fc and fc'. 

Result. — First we show the results in RPA. The or- 
bital ordering temperature, T^ PA , is identified as a tem- 
perature where the susceptibility diverges due to the con- 
dition of 



Det 1-Ux°(q,0) 



0. 
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The obtained ordering temperature monotonically in- 
creases with U . The momentum q, in which the suscepti- 
bility diverges at the highest temperature, is (71-, 7r)(= Q) 
for all values of U. 

To examine the fluctuation just above T < ? PA , the con- 
tour maps of the (zz) components of the susceptibilities 
X^ z (q,iui n = 0) for several t/'s are presented in fig. [3J 
The temperature is chosen to be T/t = Tf ph /t + 0.01 
for each value of U. In the case of small U, X^f A shows 
large intensity along (n, q z ). On the other side, in U — 5, a 
large fluctuation emerges along (q x , n). In between, a spot 
is seen around (tt,tt). That is, the anisotropy in x^F A 
is changed with the Coulomb interaction. In the RPA 
scheme, the susceptibility is explicitly given in eq. (|19p 
where x° is independent of U. Therefore, the observed 
characteristic anisotropy in x^f A is reflected from the q 
dependence of x° around T^ PA , which depends on U . 

In order to understand the origin of the above results, 
we focus on the bare susceptibility x° z (q,iw„ = 0) at the 
two representative momenta q = (tv,0) and (0, 7r). In the 
equation for X° zz {q^ * w n)j 



1 



x zz (q,iu n ) = i^2 



dk 



.z. (2tt) 



2 Cr +q , k K\ qtk (i^ n ), (25) 



where the coefficient C^ k = A\ +q A\ + - 
2rjr]' B k+q B k and a symbol 77 = ± for 77 = =p. At the 
momenta q = (tt,0) and (0,7r), this coefficient does not 
depend on k, and satisfies the relations 



C T+ Q M = <W at 9=(T|0), 



(26) 



and 



<?fe+ g>fe = l-<W at <? = (0,tt). (27) 



By using the above relations, we have the equations 

1 f rlk 

xUQ = M),0) = o / 7^ht+J0), (28) 



and 



2 7b.z. (2^) 2 k+qM 



X°,(g = (0,7r),0) = - / t^K-JO), (29) 



2 7b.z. (2tt) 2 k+qMK 

where we introduce the Lindhard function h r j^ +q k (iu n ) = 

[f{n +q )-f{El)]/{iuJ n -El +q +Ei), the Fermi distribu- 
tion function /(e) = l/[cxp/3(e — /1) + 1]. These relations 
imply that the intra-band [inter-band] nestings between 
the vertical Fermi surfaces along (— tt/2, k z ) and (7r/2, k z ) 
are concerned in x° z (<? = (tt, 0), 0) [x° zz (q = (0, tt), 0)] (see 
fig. [T]), and govern the observed anisotropy in the suscep- 
tibility. From now on, x^ z (Q = (^O^O) and X zz (q = 
(0,tt),0) are termed x°„ t ra and X° n ter> respectively. 

The temperature dependences of x? n tra an( i Xhiter are 
shown in fig. [4| In low temperatures, Xhitra shows a log- 
arithmic divergence, and Xinter i s almost constant. There 




Fig. 4: Temperature dependences of the susceptibilities x?m 
and xfntcr- 



is a crossing point, and Xinter becomes larger than Xintra 
in high temperatures. 

Schematic intra-band nestings in xpntra are shown in 
figs. [5ja) and (b). The two points, e k and e k + q , con- 
nected by q = (tt, 0) are located in the left and right 
branches in the same band in the figure. Band curva- 
tures near the Fermi surface arc opposite with each other. 
In low temperatures [see Fig. EJa)] , there are a number 
of pairs for (e k ,£k+ q ) which contributes to x? n tra with a 
fixed momentum q = (71", 0). With increasing temperature 
[see Fig. [5jb)], thermal broadening of the Fermi surfaces 
reduces the susceptibility, as well known in the nesting in 
a single-band model. 

Inter-band nestings in Xinter are a ^ so shown in figs. [5Jc) 
and (d). Band curvatures, where e k and E k + q are located, 
are parallel, in contrast to the case of Xhitra- ^ s a resv dt, in 
low temperatures, a number of pairs for (e k , e k + q ), which 
satisfy the nesting condition at q = (0,7r), are limited 
to vicinity of the Fermi surface. This is the reason why 
Xinter ^ s smaller than Xintra m l° w temperatures. With 
increasing temperature, thermal broadening of the Fermi 
surfaces reduces the susceptibilities in both the two cases 
in Xinter and Xintra- However, a reduction rate in x?„ t or 
is smaller than that in Xintra' since the nesting conditions 
in Xinter are satisfied, even when the two points, e k and 
£fc+ q , are away from the Fermi surface within an energy 
range of the order of temperature. 

Now we explain the mechanism of the characteristic U 
dependence of the anisotropy in x?z PA based on the tem- 
perature dependences of Xinter and Xintra- 111 small U, 
i.e. small Tj lPA , the dominant fluctuations are caused 
by the intra-band nestings. The susceptibility at (7r, 0) is 
larger than that at (0, 7r), because the coefficient k 
at q = (tt, 0) is finite only for the intra-band fluctuation, 
i.e. 77 = r/, as shown in eqs. (|26|) and (|27|) . With increas- 
ing the ordering temperature by increasing U, reduction of 
the intra-band contribution is more remarkable than that 
of the inter-band one, and the the susceptibility at (0,7r) 
becomes larger than that at (tt, 0). We conclude that the 
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Fig. 5: Schematic views for the band nestings, (a) [(b)] is for 
Xintra in l° w [high] temperature, and (c) [(d)] is for Xinter in 
low [high] temperature. 



observed characteristic U dependence of the anisotropy 
in x^ z PA is a consequence of i) the change of the dom- 
inant orbital fluctuation from the intra-band nesting to 
the inter-band one with increasing T, and ii) the band- 
and momentum-dependent coefficient C^_ q k . 

Next we introduce the numerical results obtained by the 
SC-RPA method. The U dependences of the susceptibility 
Xz? _RPA are shown in figs. [5] Unlike the results in fig. 
temperature is fixed at T = O.li. In the numerical cal- 
culations, a number of meshes in the first Brillouin zone 
is chosen to be 128x128, and that for the imaginary time 
between and /3 is chosen to be 1024. To calculate the 
excitation spectra, the analytic continuation is adopted as 
ie n — >• e + i5 with a small constant S = O.Olt. At T = O.li, 
the numerical iterative calculations are converged in the 
region of < U < 3.25. This implies that, at U = 3.25, 
the orbital ordering temperature is less than T = O.li 
which is lower than that in RPA. It is shown in fig. [5] that 
the U dependences of the anisotropy in xf? _RPA are s™ - 
ilar with the results in RPA (see fig. EI; with increasing 
U, large fluctuations along (jr, q z ) is changed into the ones 
along (q x ,ir). 

The results are interpreted by the one-particle excita- 
tion spectrum. This is defined by A(k,e) = V, A 1 {k,e) 
with^4 7 (fc,e) = — (l/7r)Im£ 77 (fc, ie n — s- e + iS). The result 
of the one-particle spectrum at U/t = 2.0 is presented in 
fig-Efb), as well as the density of state (DOS) in fig. EJa). 
It is seen in the inset that the band curvature is strongly 
reduced from that at U = 0. The mass enhancement is 
examined from the renormalization factor at the Fermi 
surface defined by 



z k = ^Tr 



d 

1 - — ReE(k,e) 
oe 



e=0 



(30) 



where values in the two bands are averaged. We show in 
fig. [8] the U dependence of z k at k = (tt/2,tt/2). The 
factor Zk monotonically decreases with U. The results 
imply that the effective temperature for the low energy 
fluctuation is enhanced with increasing of U as a result 
of reduction in the energy scale at vicinity of the Fermi 
level. Therefore, the interpretations for the RPA results 




Fig. 6: Contour plots of the (zz) components of the orbital sus- 
ceptibilities as a unit of t^ 1 obtained by the SC-RPA method. 
Temperatures T/t — 0.1 in all figures. Parameter values of U /t 
are chosen to be (a)0.05, (b)2, (c)3, and (d)3.25. 
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{K,Kl2) (31,0) (0,7C) (7t/2,7T) 
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Fig. 7: (a) DOS. (b) a contour map of the one-particle ex- 
citation spectrum A(k,s) as a unit of t~ 1 in the k-e plane. 
Parameter value of U/t is chosen to be 2 and T/t = 0.1. The 
dotted line in (a) and that in (b) are for DOS and the band 
dispersions, respectively, in the case of U = 0. Inset in (b) is 
an expansion at vicinity of (tt/2,tt/2). 




Fig. 8: The renormalization factor at the momentum k — 
(n/2, n/2) and a number of the states for each orbital near the 
Fermi surface defined in eq. I]31fl. 
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introduced previously are also applicable to the present 
SC-RPA results where the temperature is fixed. 

Discussion. — In this section, the numerical results 
introduced in the previous section are further discussed 
from the view point of the anisotropic transfer integral 
tj 1 in eqs. ([5]) and ([6]). Let us focus on the (zz) com- 
ponent of the bare orbital susceptibility, ~x^ zz {q, 0)j a ^ t ne 
two representative momenta q = (7r, 0) and (0,7r). As ex- 
plained in eqs. (f^5|) and (|27p . it is sufficient to consider 
the nestings between the vertical Fermi surfaces, i.e. the 
lines along (ir/2, k z ) and (— 7r/2, k z ) in fig. [ljb). Along 
these lines, the bare energy bands in eq. ([5]) are given as 
£(7r/2,fc z ) = £ (-ir/2,k z ) = + ^ Z ) cos k z , which implies 
that the a and b orbitals are not hybridized with each 
other, and the Fermi surfaces, i.e. e^/^y = £(_ n /2,k z ) = 
0, consist of the fiat &-orbital band except for the points 
(±7r/2, ±7r/2) and (±7r/2, =F7r/2). This is a consequence 
of the anisotropic transfer integral tip ; electrons in the 
b orbital do not hop along the z direction. Therefore, in 
the sense of the non-interacting electrons, the nestings be- 
tween the two vertical Fermi surfaces along (— tt/2, k z ) and 
(•7r/2, k z ), connected by the momentum q = (-7T, 0), are re- 
flected from the one-dimensional character in the 6-orbital 
electrons. 

Now we consider the effect of the Coulomb interaction U 
on the nestings. This interaction provides the hybridiza- 
tion between the a and b orbitals even on the vertical Fermi 
surfaces because of the off-diagonal matrix elements e.g. 
— U {c\c a ) c^Cb in the mean-field sense. This hybridization 
induces finite electron motions in the b orbital along the z 
direction. To check this scenario, we calculate a number of 
the states for each orbital near the Fermi surface defined 

<31) 

where e c is the cut-off energy of the order of tempera- 
ture. The numerical results with e c /t ~ 0.1 are pre- 
sented in fig. [SJ The interaction reduces more strongly 
Ab rather than A a ; the &-orbital character near the Fermi 
surface is weaken by U. This hybridization suppresses the 
one-dimensional character in the vertical Fermi surfaces, 
and the intra-band nestings which govern the fluctuation 
around (n, 0) are weaken. In other word, a ratio of the 
inter-band nesting, which has a dominant contribution at 
(0,7r), to the intra-band one is enhanced by U. 

Finally, the observed anisotropy in the orbital fluctua- 
tions is discussed with the connection to the original com- 
pass model in cq. (TTJ). It is clear from the Hamiltonian in 
eq. (TTJ) that the correlation between the z component of 
the pseudo-spins, T z , is stronger along the z direction than 
that in the x direction. This is caused by the anisotropic 
exchange interactions originating from the bond-depend 
transfer integrals in eqs. (J5j> and © in the perturbational 
sense. For example, the exchange interaction along the z 
direction is attributed to the virtual electron hopping in 
the a orbital along this direction. This fact might corre- 



spond to the present observations that Xzz(q,0) has large 
fluctuation along (q x , ir) in the case of large U where the 
contributions of the a orbital electrons play some crucial 
role as explained previously. However, in order to examine 
how the present results in the large U case are connected to 
the original compass model, further examinations, which 
are applicable to the intermediate and strong coupling re- 
gions, are required. 

In summary, we study the two-dimensional orbital- 
compass model as an itinerant electronic model. The 
Hubbard-type Hamiltonian from which the orbital com- 
pass model is reproduced is derived. This Hamiltonian 
is analyzed by the RPA and SC-RPA methods. Rod-like 
anisotropic fluctuations along {-K,q z ) in small U is changed 
into the fluctuations along (q x ,Tr) with increasing U. This 
result originates from the fact that the dominant contri- 
butions to the orbital fluctuations are changed from the 
intra-band nesting to the inter-band one with increasing 
the interaction. The present study opens a new approach 
for the orbital compass model. 
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